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Biogeochemistry is the study of how living systems in 
combination with abiotic reactions process and cycle mass 
and energy on local, regional, and global scales 
(Schlesinger, 1997). Understanding how these biogeo- 
chemical cycles function and respond to perturbations has 
become increasingly important, as anthropogenic impacts 
have significantly altered many of these cycles (Galloway 
and Cowling, 2002; Houghton et al., 2002). Blogeochemis- 
try is strongly governed by microbial processes, and it 
appears to closely follow thermodynamic constraints in that 
electron acceptor (O>, NO;, SOJ, etc.) utilization closely 
follows a priori expectations based on energetics (Vallino et 
al., 1996; Hoehler et al., 1998; Jakobsen and Postma, 1999; 
Amend and Shock, 2001). Consortiums of microorganisms 
seem to have evolved to exploit chemical potentials wher- 
ever they exist in the environment, as manifested by the 
recent discovery of anaerobic methane oxidation by sulfate 
(Boetius et al., 2000) or sulfide oxidation by nitrate (Schulz 
et al.. 1999). Three and a half billion years of natural 
selection have produced living systems capable of degrad- 
ing most Chemical potentials. We may therefore ask: If all 
ecosystem niche space is filled, is the biogeochemistry we 
observe in the environment dependent on the organisms that 
occupy that environment, or is the biogeochemistry deter- 
mined by fundamental forces, with the evolution of living 
systems being the implementation of those forces? Recent 
developments in nonequilibrium thermodynamics (NET) are 
beginning to support the latter alternative, and advances in 
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genomics are allowing us to explore microbial consortiums 
in detail. Taking advantage of ideas being suggested by 
NET, we have developed a modeling framework that views 
microbial consortiums as an inter-spectes distributed met- 
abolic network. When combined with experimental obser- 
vations, the model should help us test hypotheses that gov- 
ern how living systems function. 

The main challenge to understanding microbial biogeo- 
chemistry is understanding the complex, but mostly coop- 
erative, metabolism that develops among organisms that 
orchestrate biogeochemical cycles. Consider, for example, 
the metabolism found in the rumen of ruminant organisms 
such as cows (Madigan et al., 2000). The ecosystem of the 
rumen develops dozens of functional groups consisting of 
hundreds of species that degrade cellulose and starch to 
many intermediate organic acids and alcohols, as well as to 
CO,, methane, and hydrogen. Many of these organisms 
cannot survive without the presence of others. For instance, 
ethanol fermentation to acetate and H, is unfavorable due to 
the accumulation of H», which is also toxic to many organ- 
isms. However, if the fermenting organisms are coupled 
with methanogens (i.e., syntrophy), the overall reactions can 
proceed, and are very efficient (Jackson and McInerney, 
2002). Furthermore, these systems are controllable, hence 
predictable, as the host organism is able to utilize a food 
source that it is not capable of digesting without the organ- 
isms in its rumen. What governs the development of this 
biochemistry? Interestingly, no one organism conducts all 
these biochemical transformations. Instead, the overall me- 
tabolism of the system is distributed across hundreds of 
different microbial species. Yet the system is well coordi- 
nated due to multiple levels of organization and multiple 
levels of proliferation, not by Darwinian selection (Caldwell 
et al., 1997). Can this enigmatic coordination be explained 
by nonequilibrium thermodynamics? 
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Nonequilibrium Thermodynamics 


The application of thermodynamics to living systems 
dates back to Schrödinger (1944) and his examination of the 
creation of order from disorder. Although at first it appears 
that living systems violate the second law of thermodynam- 
ics in that they synthesize order from disorder, Schrödinger 
solved this problem by turning to nonequilibrium thermo- 
dynamics (NET). In an open system. energy flux from 
outside the system can reduce the system's internal entropy, 
which is offset by an equal or greater increase in entropy 
outside the system. Since Schrödinger, many investigators 
have extended the application of NET to living systems 
(Prigogine. 1955; Margalef, 1968; Morowitz. 1968) with 
many recent advances (Allen. 1985; Johnson. 1988: Schnei- 
der, 1988; Wiley, 1988; Choi et al., 1999; Jorgensen et al., 
2000; Toussaint and Schneider, 1988). Schneider and Kay 
(1994) succinctly describe the current restatement of the 
Second Law: 


Thermodynamic systems exhibiting temperature, pressure 
and chemical equilibrium resist movement away from their 
equilibrium states. When moved away from a local equilib- 
rium state a system will behave in a way which opposes the 
applied gradients and moves it back to its local equilibrium 
attractor. The stronger the applied gradient, the greater the 
effect of the equilibrium attractor on the system. The more a 
system is moved from equilibrium, the more sophisticated its 
mechanisms for resisting being moved from equilibrium. If 
dynamic and/or kinetic conditions permit, self-organization 
processes are to be expected. This behaviour is not sensible 
from a classical second law perspective, but is what is ex- 
pected given the restated second law. No longer is the emer- 
gence of coherent self-organizing structures a surprise, but 
rather it is an expected response of a system as it attempts to 
resist and dissipate externally applied gradients which would 
move the system away from equilibrium. 


Dissipative systems need not always operate at an opti- 
mum (Ulanowicz, 1997), because perturbations can reduce 
system organization, which will result in slower gradient 
degradation. For example. a perturbation of significant mag- 
nitude can destabilize the rumen microbial system and kill 
the host organism. lf the feed changes abruptly from cellu- 
lose (grasses) to starch (grains), the organization of the 
microbial system can collapse due to a rapid increase in the 
bacteria that produce lactic acid, thus causing acidosis. 
Interestingly, if the feed is gradually changed, destabiliza- 
tion does not occur. Such phenomena are the essence of 
self-organization, which occurs in many autocatalytic sys- 
tems (Ulanowicz, 1997). 

Simple examples of the restated second law abound. For 
example, a chamber of gas isolated from any energy gradi- 
ents soon reaches equilibrium in which the gas molecules 
are uniformly distributed (7.e., its highest entropy state). If 
the chamber is placed within a thermal gradient, the molec- 


ular gas distribution is no longer uniform, but shows in- 
creased order since a higher density of molecules will reside 
near the cold end of the chamber. Hence, the thermal 
gradient has lowered the internal entropy of the chamber. 
Once more. if the thermal gradient is sufficiently strong, a 
density-driven circulation will develop, further increasing 
system order. while reducing system entropy. This self- 
organization can also be seen in hurricanes, which are the 
organized structures that develop to facilitate degradation of 
the thermal gradient that has built up between the atmo- 
sphere and ocean over the summer, In an analogous manner, 
living systems are the organized structures that develop to 
degrade incoming solar radiation and chemical potential. 
Indeed, by examining emitted thermal radiation, Schneider 
and Kay (1994) have found that mature forests degrade 
incoming solar radiation more effectively than early succes- 
sional forests or arid land. 

The observation that ecosystem function may be gov- 
erned by laws that transcend Darwinian selection is not new. 
Many theories related to, or derived from, NET have been 
developed to describe objective functions living systems 
tend to follow; however, the choice of an appropriate ob- 
jective function still remains controversial. This is under- 
standable. since current theories of NET apply only to linear 
approximations in the neighborhood of an equilibrium (On- 
sager, 1931; Prigogine, 1978). Theories for systems far from 
equilibrium are still under development. As a result, objec- 
tive functions that ecosystems might track are numerous, 
and include optimizations of exergy (Jorgensen, 1994; 
Nielsen, 1995, Jorgensen et al., 2000) emergy (Odum, 
1983), ascendancy (Ulanowicz, 1986), power (Odum and 
Pinkerton, 1955; Odum, 1971). biomass to maintenance 
(Margalef, 1968). thermodynamic efficiency (Nielsen and 
Ulanowicz, 2000), or entropy (Prigogine and Nicolis, 1971). 
As NET theories develop, perhaps many of these observa- 
tions and theories will be collapsed or falsified. 

Our objective is to develop a modeling framework that 
can be used to test these various objective hypotheses, but 
we must first have a basis for the model. Johnson (1988) 
argues that the rotational pattern of the earth converts the 
rectilinear energy output of the sun into a pulsed energy 
input to the earth, which induces cyclical energy flows, 
resonance, and time delays. This allows the energy of the 
system to be pumped up before it decays back towards 
ground state. Metabolically, this is what we see. Autotrophs 
use incoming solar radiation to create chemical potential in 
the form of redox gradients. Material at equilibrium, H,O 
and CQ,, is pumped into a high-energy state by its conver- 
sion to O, (oxidizing) and glucose (reducing). Heterotrophs 
return this redox gradient to ground state in a cyclical 
manner. Real ecosystem biogeochemistry is more complex 
when we include abiotic reactions, anaerobic environments, 
alternate electron acceptors and donors, nutrient constraints. 
and transport limitations. Nevertheless, it is the buildup and 
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Table | 


Reactions used to model planktonic microbial community. 





Reaction 





C-fixation 
NH,” uptake 
NO, uptake 


Stoichiometry 


a CO, + 20 eee — C.H,50, + 6150 
NH,” + 0.58 C,H,.0, + 0.034 SO,7" + 0.9e © AA + 1.9 H,O 
NO, + 10H* + 0.58 C,H,,.0, + 0.034 SO, + 8.9 e7 — AA + 4.9 H.O 


Nitrification NH, ANOT WO, + |IOH +8 es 

N, fixation aN, + 4H* + 0.58 C,H,,0, + 0.034 SO + 39e° — AA + 19H,0 
Photosynthesis 2H,O>O,¢e8H tde 

Respiration O, +™4H™ Fe — 2H,0 

Protein AA + (Protein), © (Protein),,,; + HO 

Chlorophyll 4 AA + 6.8 C,H,,0, + Mg** + 46.2 Ht + 48.6 e7 < Chl + 42.3 H.O + 0.14 SO, 
Lipids 8 CH 20, + PO,” + 92H" + 89e < CisHo OP + 46 H20 

Glycogen C6H 206 + (CoH inOs)a S (CoH iuOs)ne + H20 

DNA/RNA iene PO, & DNA + 0.73 CH -0 + OM SO2 + 19HiO+ Io > 75m 





Abbreviations: AA, amino acids (C3 H7,0O,7NSy 934° °°”): 
(Cg 6H 1408N,P-). Reactions are either reversible (<>) or irreversible (—). 


decay of redox potential via a distributed metabolism that 
forms the cornerstone of our approach. We focus on micro- 
bial systems because these systems exhibit the greatest 
degree of metabolic capacity, are responsible for the major- 
ity of biogeochemistry on earth, display fast dynamics that 
allows for practical experiments, and are Jess susceptible to 
loss of diversily. 


Chi, chlorophyll) (C5.H>,0.N,Meg); Protein (C,.H.,O,,NSiee DN 


Thermodynamically Constrained Metabolic 
Biogeochemical Model 


A traditional reductionist biogeochemical model would 
include differential equations for growth of each microbial 
functional group, equations for Monod-type growth kinet- 
ics, and numerous conditional statements to direct the use of 
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%20, + 2H* + 2e° > H,O 
NO, + 6H* + 5e > “AN, + 3H,O0 
NO, + 2H* + 2e-— NO, + H,O 

SO,* + 10H* + 8e° > H,S + 4H,O0 
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CO, + 3H* + 4e° — CH,O + OH 
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Figure l. Conceptnalization of the metabolic biogeochemistry model. Half reactions lead to production of 
protein (and other building-block constituents), which is then allocated to reactions governed by an optimization 
function. Abiotic reactions are incorporated with standard kinetics. 
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Tahle 2 


Governing equations for the thermodynamically constrained 


biogeochemistry model 





Coastraral Equation 
Free energy aG, 
Redox Dre = 0 
7 y TÈR (n) 
Biomass composition GN = TR) = CN 
Reaction kinetics pea! one Py neg 
z Weer Sane i aE, dE, dE, 
c ar 9 = í —— = - = — 
nzyme allocation E r and ly Ai d 
| PAR(1) | [Chl] 
Photosynthesis ip = 7 NEA ee 
LP Problem Maximize SR EVIER e E) 
nE 


f 


Definitions: r, are reachion rates (Table 1); e, are electron-pair transfers 
associated with reaction i; ŒE, 1s conceatralion of enzyme allocated to 
reaction i, and E; is the total protein concentration: R© and R” are the 
subset of reactions (7,) that lead to C and N accumulation in living biomass, 
respectively; f(S) are functions that describe uptake kinetics of substrate S. 
such as Michaelis-Menten; PAR(z) is the photosynthetic active radiation; /: 
is water column depth; and r, is the rate of photosynthesis (Table 1). 
Solution of the linear programming (LP) problem gives r, and Æ, at time £, 
which are used in standard C and N conservation equations to obtain 
solution of stale variables over time. See Vallino et al. (1996) for a more 
thorough development. 
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Figure 2. 


the various electron donor and acceptor resources that exist 
within the environment (Koelmans et al., 2001). While this 
approach has been useful for well-defined laboratory exper- 
iments with a few species growing on a limited number of 
well-defined substrates, it is not practical for extension to 
more diverse microbial ecosystems with numerous or ill- 
defined substrates. Reductionist-based biological models 
fail to incorporate the governing laws that define living 
systems (Lawton, 1999); the models are based solely on 
empirical observations. Consequently, these models are 
brittle and often fail as the system’s state changes signifi- 
cantly over time and space (Vallino, 2000). 

To develop a robust model that can predict microbially 
governed biogeochemistry tn spatially and temporally di- 
verse environments, a more holistic, systems-based perspec- 
tive must be taken. Our governing philosophy ts that living 
systems synthesize and allocate metabolic capability in such 
a way as to optimally utihze available resources in the 
environment as governed by NET. What we seek to deter- 
mine is the nature of the objective function that living 
systems tend to follow, and what causes living systems to 
diverge from this function. 

This optimization-based approach was first developed in 
a thermodynamically constrained metabolic framework to 
examine bacterial utilization of dissolved organic matter 
(Vallino et al., 1996). However, this model still uses an 
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Metabolic model. Changes in resource concentrations of (a) ammonium aad (b) nitrate, and 


accumulation of biological structure (c) protein and (d) chlorophyll over the course of the simulation. 
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(b) NO, Uptake 







1.2 


Rxn Rate (.M/d) 
Rxn Rate (uM/d) 


(d) 


Enzyme Conc. (uM) 
Rxn Rate (uM/d) 


1 2 
Time (d) 


Reaction rates (solid line) and allocation of protein (enzyme concentration, dashed line) tor 


reactions involving (a) ammonium uptake, (b) nitrate uptake, (c) N, fixation, and (d) chlorophyll synthesis during 
the course of the simulation. Abbreviations: Gle, glucose; AA, amino acids, Chl, chlorophyll. 


organismal approach, in that the model tracks bacterial 
biomass. We increase the applicability of our modeling 
approach by removing emphasis on synthesizing bacterial 
biomass and placing it instead on synthesis of metabolic 
capability exhibited by the whole ecosystem. The model 
consists of a set of metabolic half-reactions that represents 
the major metabolic capability of a planktonic ecosystem 
(Table 1). But instead of synthesizing bacteria, reactions 
produce protein, chlorophyll, storage compounds, and other 
fundamental building materials observed in living systems 
(Fig. 1). These building materials represent those summed 
over all organisms in the ecosystem, not any one particular 
organism. Indeed, organisms are not directly modeled. 
Newly synthesized protein is then allocated to those meta- 
bolic reactions that optimize the specified objective criteria, 
while enzymes no longer in use can be degraded back into 
constituent amino acids (Fig. 1). A linear programming (LP) 
problem is used to determine the reaction rates (7;) and 
enzyme concentrations (£,) that maximize a given objective 
function, subject to fundamental constraints, such as energy, 
redox, composition, kinetics. and light-capturing capabili- 
ties (Table 2). Although the model does not distinguish 
species in a classic sense, it does from a functional perspec- 
tive. As environmental conditions change, so do allocations 
of resources to metabolic reactions. Real systems accom- 


plish this same objective via relative changes in species 
abundances and magnitude of gene expression. 

As an example of the model. we simulate a marine 
phytoplankton bloom, where metabolic reactions associated 
with ammonium (NH,°) and nitrate (NO; ) uptake, N,- 
fixation, carbon dioxide (CO,) fixation, and biomass syn- 
thesis (protein and chlorophyll) are included in the model 
(Table 1). The optimization goal chosen was maximizing 
the rate of biomass synthesis, though others could be for- 
mulated. Resources made available were 5 uM NH,*. 5 uM 
NO, , atmospheric N,, and light. The model simulation 
proceeds by preferentially consuming NH,” over NO, 
(Fig. 2a,b), which is evident by the allocation of protein (in 
the form of enzyme) to NH,” uptake (Fig. 3a). but not to 
NO, uptake nor N, fixation (Fig. 3b.c). There is also a 
strong initial allocation of protein to chlorophyll synthesis 
(Figs. 2d, 3d), but this protein is rapidly reallocated after 
0.5 d due to diminished returns on the investment in light 
harvesting capacity (/.e., chlorophyll), which saturates at 
high chlorophyll concentration (Fig. 3d). As NH, becomes 
exhausted (Fig. 2a), protein is reallocated from NH,’ 
to NO, “uptake (Fig. 3a, b). Subsequently, as NO,  be- 
comes depleted (Fig 2b), protein is allocated to N, fixation 
(Fig. 3c). 
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Conclusions 


If nonequilibrium thermodynamics governs biogeochem- 
isiry, our Melabolic modeling approach represents a more 
direct means of cupluring ecosystem dynamics than classic. 
orgamsmal-based approaches. The approach also predicts 
how whole-system genomic transcription and translation 
should proceed, which can be compared to actual systems 
using techniques currently being advanced in molecular 
biology. Because the metabolic ecosystem model is based 
on fundimental governing equations, 1t should prove more 
robust and have a grealer operaling range than organismal- 
based models. 
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